Spin-Orbit-Mediated Spin Relaxation in Graphene 
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We investigate how spins relax in intrinsic graphene. The spin-orbit coupling arises from the band 
structure and is enhanced by ripples. The orbital motion is influenced by scattering centers and 
ripple-induced gauge fields. Spin relaxation due to Elliot- Yafet and Dyakonov-Perel mechanisms 
and gauge fields in combination with spin-orbit coupling are discussed. In intrinsic graphene, the 
Dyakonov-Perel mechanism and spin flip due to gauge fields dominate and the spin-fiip relaxation 
time is inversely proportional to the elastic scattering time. The spin-relaxation anisotropy depends 
on an intricate competition between these mechanisms. Experimental consequences are discussed. 



Graphene can be useful in future advanced applications 
because of the reduced dimensionality, the long mean free 
paths and phase coherence lengths, and the control of 
the number of carriers Among possible applications, 
graphene is investigated as a material for spintronic de- 
vices [1,I3,I3,H,[1,0,S|- Spintronics aims to inject, detect, 
and manipulate the electron spin in electronic devices. 

Spin manipulation via the spin-orbit (SO) coupling has 
been extensively discussed in semiconductors and metals 

The SO coupling enables electric, and not just mag- 
netic, control of the spin In two dimensional (2D) 
semiconducting structures, inversion asymmetry results 
in the Rashba SO coupling Additionally, bulk in- 

version asymmetry in A3B5 compounds causes the Dres- 
selhaus SO coupling [12| . Device performance is limited 
by spin relaxation and understanding its origin enables 
enhanced spin control. Two mechanisms of spin relax- 
ation discussed in the literature [3] j the EUiof- Yafet 

Ti ] and Dyakonov-Perel 1^, 13] mechanisms, can be 



relevant in graphene. 

Elliof- Yafet (EY) spin relaxation is related to how 
the spi n changes its direction during a scattering event 



l4 This is possible because the SO coupling pro- 



duces electronic wave functions that are admixtures of 
spin and orbital angular momentum. Dyakonov-Perel 
(DP) ^\n\ spin relaxation is related to spin preces- 
sion between scattering events by the effective (Zeeman) 
magnetic field induced by the SO coupling. This SO 
induced effective (Zeeman) magnetic field changes direc- 
tion during scattering. In the EY mechanism, the spin 
relaxation time is proportional to the elastic scattering 
time TeijT^^ oc Tei, whereas the dependence is opposite 
oc (toi)""'^ for the DP mechanism. This qualitative 
difference allows detection of these two competing mech- 
anisms in disordered samples. 

Recently, spin transport and spin relaxation were stud- 
ied in relatively dirty graphene samples [3, [13] . A spin 
relaxation length Agf ^ 2/im was measured at room tem- 
perature and it was indicated that Agf is proportional to 
the elastic mean free pathLi, suggesting the EY mech- 
anism to be dominant [3, [l3] . The measured spin relax- 
ation length is weakly anisotropic, such that spins "out 




FIG. 1: (Color online). (a) Dyakonov-Perel spin relaxation: 
The SO coupling induces a momentum dependent effective 
field _B|| which changes direction randomly after scattering 
events, leading to spin relaxation. B\\ is "in plane" so spins 
directed perpendicular to the plane relax faster than as spins 
in the plane, (b) Gauge- Field spin relaxation: Gauge fields 
due to ripples (light and gray blue areas) induce an effective 
field B±, which for a finite SO coupling leads to spin relax- 
ation. Bi_ is "out of plane" so spins directed perpendicular 
to the plane relax slower than spins in the plane. 



of plane" relax 20% faster than spins "in plane" [18j . 
These experiments motivate a study to see if known spin- 
relaxation mechanisms, or possibly novel effects, domi- 
nate spin scattering in graphene. 

In this Letter, we consider spin relaxation in intrin- 
sic graphene arising from three ingredients: (i) The SO 
interaction arises from the band structure and can be 
enhanced by graphene corrugations, (ii) spin isotropic 
scattering centers cause momentum relaxation, (iii) topo- 
logical lattice disorder induce gauge fields that change 
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the orbital motion. First, we study the Elliot- Yafet and 
Dyakonov-Perel mechanisms for this model ignoring ef- 
fects of gauge fields. Second, we find that a unique in- 
terplay of SO and gauge fields (GF) due to topological 
disorder causes spin relaxation. Our main findings are 
that DP and GF mechanisms are comparable and domi- 
nate the EY mechanism in intrinsic graphene. Interest- 
ingly, the DP mechanism implies that spins out of plane 
relax twice as fast as spins in plane, but GF exhibits the 
opposite behavior, spins in plane relax faster than spins 
out of plane (see Fig. 1 for details) . The spin-relaxation 
anisotropy depends on the intricate competition between 
the DP mechanism and GF. Our results are valid for rel- 
atively clean graphene samples where e.g. adatoms do 
not alter the SO coupling. 

Spin- orbit coupling.- There are Rashba and Dressel- 
haus "like" SO interactions in graphene and we disre- 
gard the latter , [l^ [iSl • The total Hamiltonian reads 



H = ztvpk {fi ■ a) — — ((T X s) 



(1) 



where = ) corresponds to the K{K') point, the 
spinor basis for K is '^'k = (^f ^i' ^l)^ ^^'^ ^' 
the components are reversed, 5'^' = {B^ , Bi, A-^ , Ai)'^ , 
A and B denote the lattice sites and t and I the electron 
spin. With this choice, the SO coupling is identical for K 

and K'. n = (cos 6*, sin 6*), tan6' = ky/k^, k = ^k'^ + ky 
is quasi-particle momentum with respect to the K{K') 



corners of the hexagonal Brillouin zone, A = Ac 



-A, 



[lit is the dominant Rashba SO coupling constant in- 
duced by curvature and/or an external electric fields [l9j 
and a = {ax, (Jy), s — {sx, Sy) correspond to Pauli matri- 
ces in sublattice and spin space, respectively. The eigen- 
states of (dl for the K valley are 



C-s 



T) 



ik-r 



with energies Ek^s = sA/2 + (,D, D ^ y/ivpk)'^ -F A2/4, 
where s = ± corresponds to the "f (|) spin states and 
^ = ± denotes the pseudospin degeneracy, Ck^s — 



y/Ek^s/v2D, Cs = Ck^s and N^s = ^/2E^s / {v^k) . 
These eigenstates are polarized in plane of the graphene 
layer. The spin precession length is Zprcc = 2-kvf/ A- 

Elliot-Yafet mechanism.- The "Rashba" SO coupling 
can change the spin orientation during a scattering event 
as is typical for the EY mechanism. We study this effect 
by a decomposition into partial waves with a well defined 
orbital angular momentum discussed in [2lj. Neglecting 
mixing of the K and K' valleys, an incoming wave with 
total angular momentum, C = Tidg ± (7^/2 Sz/2, is an 



eigenstate of Eq. ((T|): 



c_ J„(fcr)e' 

C+Jn+l 

c+J„+i(fcr)e*("+i)« 
c_J„+2(fcr)e*("+2)» 



(2) 



with c± as defined earlier and J,i{x) is a Bessel function. 
We analyze here weak scatterers, where the elastic scat- 
tering rate scales as r " j""^ ~ Ep. Strong and resonant 
Coulomb scatterers [21| induces a similar change in the 
spin orientation. We approximate the potential by a step 
function, v{r) = Vq [1 — 0(r — R)]. The wave function 
inside the potential well is a superposition of two radial 
waves, finite at the origin and with different spin orienta- 
tions. In graphene, Rashba spin coupling entangles spin 
and pseudospin and the change of spin in a scattering 
event depends on the evolution of the pseudospin. Un- 
like conventional semiconductors f9], the spins in differ- 
ent angular momentum channels differ, complicating the 
definition of the amount of spin relaxation in a scatter- 
ing event. From the given incoming wave 5'*" with spin 
parallel to the momentum k and incident angle 9, there 



are two possible outgoing waves {r, - 



which satisfy conservation of energy and momentum. 
These can be written in a similar way as Eq.([2), where 
0' ~ -0 + Acot{e)/{vFk). We then define 



' On) 



(3) 



where r„(ro„) and r^{rQ^) are the scattering ampli- 
tudes for a given angular momentum channel n with 
(without) SO coupling and for the two possible outgo- 
ing waves vl(°"*^ \[/o"* respectively. SO coupling changes 
the wavevector for one of the reflected waves, k' « 
k — A/vf . The leading reflection coefficient, rn=o, in 
the absence of SO coupling depends on wavevector as 
ro{k) - VokR^/vF^, so that r„=o(A:') - rn=o{k) ~ 
VqAR'^ / {vf)'^ and 5 ~ A/(Di?fc). S vanishes if the spin 
is conserved in the scattering event. If the changes in- 
duced by a finite A are small, this quantity should be 
proportional to the change in spin orientation during 
the scattering process. The change in spin orientation 
at each collision is A/ (vFkF). The total change of the 
spin after Ncoii colhsions is of order \/Nco\\A{vFkF)~^ ■ 
Dephasing occurs when ^ Nco\iA{v pkp)^^ = 2tt after a 
time Tso = TciA'coii, where Tei is the elastic scattering 
time, and the Elliot- Yafot spin relaxation time is ~ 
{vFkFf/A^ X (Ici/vf) ^ {vFkpf/A'^ x t^i. The spin 
diffusion length is related to the spin relaxation time by 
Aso = VDtso, where D = v%t/2 so A^^ ~ le\{vFkF)/A. 

The D'yakonov-Perel' mechanism. - Between scattering 
events, the Rashba SO coupling acts as an effective mag- 
netic field in the plane, B^\{a) = A{a x ez)/2, on the 
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spins 1^, 17 1 . The spin dynamics is a result of spins pre- 
cessing in a fluctuating in-plane magnetic field governed 
by the momentum. Elastic scattering randomizes the 
momentum and the associated magnetic field. Averag- 
ing over many collisions, the spin orientation becomes 
random after a time f9] 



(4) 



e.g. the spin- relaxation time is inversely proportional to 
the elastic scattering time. The DP spin-diffusion length 
is independent of the mean free path A^-^" = vfI{V21^) = 
^proc/ (2\/27r). 

Effective gauge field.- Topological lattice defects, 
strains, and curvature change the hopping integrals be- 
tween the lattice sites. These effects are captured by an 
induced effective gauge field A{r) which deflect electrons 
and change the electronic states at low energies P, |2§| . A 
is related to the strain tensor Uij (f) describing topologi- 
cal lattice disorder , Ax ~ [^{f^xx — 'U'{r)yy\, Ay ~ u{r)xy 
and give rise to a random out-of-plane magnetic field 
B±_ = [Vx^(r)]2 We will demonstrate that this out- 
of-plane orbital magnetic field together with the Rashba 
SO coupling can induce an out-of-plane spin polarization. 

Let us first present a semiclassical argument for why 
the spin polarization changes from in plane to out of 
plane as the gauge field increases. S_l causes the electrons 
to move in cyclotron orbits with a radius — 2e/ {evB^^^ 
where v is the velocity and e is the energy. A qualita- 
tive change in the electronic states occurs when the cy- 
clotron orbit is smaller than the spin precession length, 
e.g. when 27rrc — Iprcc- Using v = vp, the graphene dis- 
persion e = vpkp, and introducing the magnetic length 
Ib = \/s/l3±, there is a transition when the magnetic 
length is shorter than 



(5) 



where Xp — 2t: /kp is the Fermi wavelength. In the low 
field regime, Ib 3> in-plane momentum is a good 
quantum number and via the SO interactions spins are 
polarized in plane. In contrast, in the strong field regime, 
Ib ^ 1%^ out-of-plane angular momentum is a good quan- 
tum number, and the SO coupling changes the polariza- 
tion of the states to out of plane. 

We will now carry out a quantum mechanical calcu- 
lation which will give more details and confirm our ar- 
guments above with the threshold value Eq. ([5]) for the 
gauge field. Inserting the gauge fields in the Hamiltonian 
in Eq. ([1]), results in a 4 x 4 Hamiltonian for one valley 



/ vpii \ 

-iA vpil 

upnt iA 

VF& / 



where II = Px — iPy, P = —iV + eA and a similar 
Hamiltonian can be written for K' , where the sign of 
the gauge field is reversed. First, we consider the so- 
lution of the Dirac equation when A = 0. The results 
are gauge-invariant and we choose the Landau gauge 
Ax = 0,Ay = B±x and we consider a homogenous mag- 
netic field to illustrate the main effect. The wave func- 
tions are 



Ti<l>n~i(x - Xo) 

cj)n{x - Xo) 



T,i), 



(7) 



with eigenenergy e„ — ±VF\/2\n\/lB in terms of solutions 
of a particle in an harmonic oscillator potential — 
xq). The two components of *&n,T,j, correspond to the 
amplitudes in the two graphene sublattices and xo = klg 
is the Landau level guiding center. 

The Rashba SO coupling induces an interaction be- 
tween electrons with spin up in one sublattice and elec- 
trons with spin down in the other sublattice. There- 
fore it is convenient to express the Hamiltonian ^ in 
the basis <^„_i[l,0f I T), </'n[0,l]^| T), 0n[l,O]^| i). 



n 



[0,lf I i),n^HvFV2/lB: 



iy/lnl 


iA 






-iA 




V2|" + l| 
-iy^\n+l\ 



(8) 



where A = AIb / {vpV2). The eigenenergies of Ti are: 



= cr[l + A^ + 2n + sW 1 + A2 +4nA2]/2, 



where, a = + {a = —) denotes electron (hole) like excita- 
tions and s = + (s = —) denotes spin. Let us consider the 
expectation value of the out-of-plane spin polarization of 
these states, pz. For ti = there are three physical states 
of which one has polarization (fi/2)(l — A'^){1 + A^) and 
two have polarizations — (ft/2)/(l -I- A^). For all states 
where n > I the polarizations are 



A2) +4nA2]-i/2 



(9) 



The spin polarization differs between the lowest and high- 
est Landau levels and the transition roughly occurs when 
4y^A « 1. Transport is governed by states at the Fermi 
energy e„ « vp ^J2\n\/lB vpkp and the condition 
Ay^A w 1 can be rewritten in terms of a critical value 
for the magnetic length, which exactly agrees with our 
semiclassical estimate for Ig in Eq. ([5]). Equation Q 
demonstrates that when Ib ^ Ib^ the states are fully 
out-of-plane polarized, but the out-of-plane polarization 

(c) 

vanishes when the gauge field is weak Ib 3> Ib ■ The en- 
(6) ergy splitting is of the order A for weak gauge fields and 



reduced by a factor 



for stronger gauge fields. 
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This change in the polarization direction of the eigen- 
states with increasing gauge fields has consequences for 
the spin relaxation. When the gauge fields vanish, the 
effective magnetic field is in plane so that spins out of 
plane relax twice as fast as spins in plane. In the regime 

(c) 

around Ib — Ig , the in-plane and out-of-plane compo- 
nents of the effective magnetic field are comparable and 
we expect the spin-relaxation anisotropy to be reduced 
and eventually exhibit the opposite behavior, spins in 
plane relax faster than spins out of plane. 

Experimental consequences.- We use typical parame- 
ters for graphene, hvp = 5.3 x eV m and the en- 
hanced Rashba coupling for a ripple of radius R ~ lOOnm 
is Ai{=ioonm = 1-7 X lQ~^eV [l9|. The Fermi wavelength 
depends on the electron doping n, Xp = y^Tr/n . In 
Rcf. [3] n « 3.6 X lO^^m^^ so Af « 13nm and ^ci = 36 
nm, but considerably larger mean free paths have been 
measured and should be expected in clean systems in the 
future. 



Comparing the DP and EY relaxation mechanisms, 
~ [v^kpf/^l^T, and r,^^ « WA^^, we find 
{kploif. Typically (/cf^ci)^ > 
i 300, so the DP mechanism is 



^DP 



' so / so 

1, e.g in Ref. 
much more important that the EY mechanism. 

The experimental trend that the spin-relaxation length 
is proportional to the mean free path [l8| is encourag- 
ing since it suggests clean systems should have a very long 
spin-relaxation length. However, this is at odds with our 
results for intrinsic graphene where the spin-relaxation 
length only weakly depends on the mean free path. The 
good news is that the computed spin relaxation length 
is long, we find Iprec = '^'^^vp / ^R=wonm = 190^m 
and hence the DP spin relaxation length is X^^^ — 
Lrec/(27r\/2) ^ 20/im. On the other hand, in Ref. 
[3| As/ ^ 1.3 — 2/xm. Our theory neither quantitively 
nor qualitatively explains the experiments in Ref. [3|. 
Adatoms such as2 hydrogen that locally enhance the SO 
interaction could be responsible for the discrepancy since 
the mobility of the experimental samples is relatively 

lowfii. 

Our theory applies to cleaner, intrinsic graphene, pos- 
sibly with less adatoms, where we predict novel spin- 
relaxation anisotropy effects. We expect it is possible 
to obtain longer spin-relaxation lengths than in Ref. 
01 t)ut that one eventually will enter the intrinsic 
regime where the DP mechanism prevents a further en- 
hancement. Furthermore, in intrinsic graphene, gauge 
fields due to ripples are important. From the parameters 
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above, we find a threshold magnetic length of I 
nm. Surface corrugations give rise to effective magnetic 
lengths of the order Ib ~ lOOnm 221 so already at these 
electron densities gauge fields are important and reduce 
the spin relaxation with respect to DP mechanism. Lower 
electron densities, which are feasible, should enhance the 
effects of the gauge fields. In ultra clean systems, /ci > ^s, 



we expect t^J' vfI~b I ^ and the spin-relaxation time 
saturates. We encourage experiments in clean graphene, 
which will further elucidate the interplay between SO 
coupling, momentum relaxation, and gauge fields. 
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